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The phase diagram, (T, p), of a finite, constrained, and 
classical system is built from the analysis of cluster distribu- 
tions in phase and configurational space. The obtained phase 
diagram can be split in three regions. One, low density limit, 
in which first order phase transition features can be observed. 
Another one, corresponding to the high density regime, in 
which fragments in phase space display critical behavior of 
3D-Ising universality class type. And an intermediate den- 
sity region, in which power-laws are displayed but can not be 
associated to the abovementioned universality class. 

PACS number(s): 25.70 -z, 25.70.Mn, 25.70.Pq, 
02.70.Ns 



I. INTRODUCTION 

The multifragmentation phenomenon that takes place 
in nuclei with excitation energies above 2 MeV/nucleon 
has been one of the central issues of the nuclear commu- 
nity during the past two decades. Many experimental 
and theoretical efforts have been concentrated towards 
the understanding of the mechanisms involved in that 
process. In particular, one feature that triggered the in- 
terest on the field was the important detected production 
of intermediate mass fragments (IMF's). Highly excited 
nuclei break up in many IMF's in Fermi energy range 
reactions. This feature, along with the fact that early 
calculations of caloric curves showed an approximately 
constant behavior, was interpreted by many groups as a 
signature of a phase transition taking place in finite nu- 
clear systems 0-H]. In recent contributions the behavior 
of the microcanonical heat capacity was used to experi- 
mentally characterize the transition as a first order one 
||. Within this picture, the latent heat can be associated 
with the transformation between a Fermi liquid and a gas 
phase, composed by light particles and free nucleons. 

Despite that the familiar liquid-gas transition frame- 
work seems to be appropriate to deal with the nuclear 
case, there are some peculiarities that are worth to be 
considered. One key point usually disregarded is that, by 
its own nature, multifragmentation in nuclear reactions 
should be a priori analyzed as a nonequilibrium process 
(see )• Even if an equilibrium scenario is adopted, a 
major point is that we are dealing with a phase transition 
occurring in a finite system. Therefore, several thermo- 
dynamical features, e. g. the entropy extensivity, cannot 
be taken for granted and the agreement between differ- 
ent statistical ensembles predictions cannot be invoked 
anymore JlO) ] . A working hypothesis, usually made by 
statistical models in order to apply a tractable global 



equilibrium picture, is that a freeze-out volume can be 
defined inside of which the existence of an equilibrated 
ensemble of clusters can be assumed. In these models, 
the behavior of thermodynamics quantities is closely re- 
lated to the way the system is partitionated into clusters 
giving rise to internal surfaces (isolated drops). 

The aim of the present paper is to analyze the relation- 
ship between cluster distributions inside a fixed volume 
and the respective thermodynamic description. In order 
to do that a confined system that interacts via a two- 
body Lennard- Jones potential is studied. A molecular 
dynamics, MD, approach is used. The analogy between 
the nuclear force and the van der Waals interaction sup- 
ports the use of this simplified classical model to obtain 
qualitatively meaningful results. 

Given the MD microcanonical description, a connec- 
tion with thermodynamical quantities can be established 
via well known principles from mechanics. For instance, 
the pressure and temperature can be estimated from the 
generalized equipartition theorem |Tl[| , whereas the spe- 
cific heat can be related to kinetic energy fluctuations 
|]l2| . On the other hand, a complementary analysis of 
the microscopic correlations at a "nucleon" level of de- 
scription can also be considered. To that end, different 
fragment-recognition algorithms can be used in order to 
unveil different particle correlation properties. 

Previous studies [jl3| have already dealt with this 
'cluster structure-thcrmodynamic description' mapping. 
They made used of a cluster definition a la Coniglio- 
Klein, and were mainly focused on the system behavior 
in the supercritical phase {p > p Cl T > T c ). 

In this contribution we consider two alternative frag- 
ment definitions that, differently from Coniglio-Klein 
clusters (that in the present paper will be called MSTE 
clusters, as will be explained in Sec. Ill), are built in 
well defined and physically meaningful spaces. The first 
one, associated with the so called minimum spanning 
tree fragment recognition method, is based on configu- 
rational information (MST clusters) H]. The second 
one uses complete phase space information in order to 
define a fragment set according to the most bound den- 
sity fluctuation in phase space []l6| (ECRA clusters). Us- 
ing this fragment characterization as a fundamental piece 
of information, consistent phase diagrams are built from 
scratch for the analyzed finite systems. In addition, the 
relationship between critical signatures in the fragment 
distributions and the corresponding system localization 
in the phase diagram is established. 

This paper is organized as follows. In Section || we will 
describe the model used in our simulations. Section III is 
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devoted to a characterization of the used cluster defini- 
tions. A description of expected signals associated with 
phase transitions in finite systems is given in Section IV . 
In Section the study of phase diagram is presented. 
Finally, in Section [vj, conclusions are drawn. 



II. THE MODEL 

The system under study is composed by excited drops 
made up of particles interacting via a 6-12 Lennard Jones 
potential, which reads: 



V(r) 
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We took the cut-off radius as r c — 3er. Energies 
and distances are measured in units of the potential 
well (e) and the distance at which the potential changes 
sign (er), respectively while the unit of time used is: 
t = y/a 2 m/4:8e. 

In order to constrain the dynamics we used a spherical 
confining 'wall'. The considered external potential be- 
haves like V wa ii ~ (r — r wa u) with a cut off distance 
r cut = Iff j where it smoothly became zero along with 
its first derivative. The set of classical equations of mo- 
tion were integrated using the well known velocity Verlet 
algorithm fr7[] , taking tint = 0.002i Q as the integration 
time step. Once the transient behavior was over a micro- 
canonical sampling of particle configurations every 5io up 
to a final time of 140000£o was performed. 



III. FRAGMENT DEFINITIONS 

The simplest and more intuitive cluster definition is 
based on correlations in configuration space: a particle i 
belongs to a cluster C if there is another particle j that 
belongs to C and |r- 
called the clusterization radius 

We set r c i = r c = 3<7. The algorithm that recognizes 
these clusters is known as "Minimum Spanning Tree" 
(MST). In this method only correlations in q-space are 
used, neglecting completely the effect of momentum. As 
was shown in Ref. fl4], MST clusters give incorrect infor- 
mation about the meaningful fragment structure of the 
system for dense configurations. However, an interesting 
point to be notice is that it can still provide useful in- 
formation about the limit imposed by the constraining 
finite volume to the formation of well defined fragments 
in configurational space. As r c ; = r cut , no inter-cluster 
interaction exists, so cluster surfaces can be univocally 
defined. 

An extension of the MST is the "Minimum Spanning 
Tree in Energy" space (MSTE) algorithm. In this case, a 
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where e# = Vfcj) + (p, — Pj) 2 /4/z, and \x is the reduced 
mass of the pair {i,j}- This cluster definition resem- 
bles the clusterization prescription adopted by Coniglio 
and Klein |L8| and was the one used by Campi et al. 
in Ref. |]|. The MSTE algorithm searches for config- 
urational correlations between particles considering the 
relative momenta of particle pairs, and typically recog- 
nizes fragments earlier than MST in highly excited un- 
constrained systems H , E5| . 

Finally, a more robust cluster definition is based on 
the system "most bound density fluctuation" (MBDF) in 
phase space [^6|. The MBDF is composed by the set of 
clusters {Ci] for which the sum, E{d}, of the fragment 
internal energies attains its minimum value: 
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Kj m is the kinetic energy of particle j measured in 
the center of mass frame of the cluster which contains 
particle j, and Vij stands for the inter-particle potential. 

The algorithm that finds the MBDF is known as 
the "Early Cluster Recognition Algorithm" (ECRA). It 
searches for simultaneously well correlated structures in 
both, q, and p space, via the minimization of the poten- 
tial and the kinetic terms respectively. 

The ECRA algorithm has been used extensively in 
many studies of free expanding fragmenting systems 
and has helped to discover that excited drops break 
very early in the evolution. In addition, in a recent contri- 
bution || it was shown that ECRA clusters are also suit- 
able to describe the fragmentation transition that takes 
place in volume constrained systems. 



IV. CHARACTERIZING THE TRANSITION 

As mentioned in the introduction, several observables 
can be studied in order to analyze phase transitions phe- 
nomena occurring in finite systems. They usually involve 
the behavior of caloric curves P,p|, ^0| J23| , specific heats 
^^J, kinetic energy fluctuations iPtpL fragment mass 
distributions P§,E6[ , and critical exponents p6[ . 



A. Thermodynamical Description 

One of the most relevant quantities in the study of 
fragmenting systems, either constrained or free to ex- 
pand is the caloric curve (CC). The CC is defined as the 
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functional relationship between the system energy and 
its temperature, given by: 



T(E) 



3(N-1) 



<K>, 



(4) 



being N the number of particles, and < K >e the mean 
kinetic energy averaged over a fixed total energy MD sim- 
ulation (9). 

In Figure |l|(a) the CC's are shown for the following 



densities: p h = 0.55a 
and p l = O.Oler- 3 . 
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FIG. 1. The caloric curve is shown in panel (a). In pan- 
els (b) and (c) the potential energy (V) and its root mean 
squared deviation as a function of the system total energy are 
displayed in panels (b) and (c) respectively. The following 
density values: p — 0.01,0.07,0.35, and 0.55<r~ 3 , displayed 
with full triangles, diamonds, squares, and circles respectively, 
were considered. 

Different behaviors can be recognized. For the more 
diluted case, p l (full triangles), the corresponding CC 
displays a loop which ends in a linearly increasing tem- 
perature line, which we refer as vapor branch. Taking 
into account that, within the microcanonical ensemble, 
the specific heat is defined as: 
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it is clear that for this case negative values of C will 
be found in the range 0.1e < E < 0.6e . A negative 
value of the derivative of the temperature as a function 
of the energy reflects an 'anomalous' behavior of the sys- 
tem entropy for the respective energy range. A convex in- 
truder in S(E), prohibited in the thermodynamical limit, 
arises as a consequence of the finiteness of the system and 
the corresponding lack of extensivity of thermodynami- 
cal quantities like S. This signal is expected in first order 
phase transitions, and is associated to a negative branch 
of the heat capacity between two poles |I(J . Taking into 
account the relationship between kinetic energy fluctua- 
tions and the system specific heat ffsU^T 



(6) 



it can be seen that negative values of the specific heat 
appear whenever (5K) 2 get larger than the canonical ex- 



pectation: 
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This behavior has already been verified 



in Ref. §. 

As the density is increased the loop is washed away and 
it is replaced by an inflection point (p" 1 , full diamonds in 
Fig.[|(a)). This corresponds to the merging of the specific 
heat negative poles into a single local maximum. Finally, 
for higher densities, p c and p h (filled diamonds, and cir- 
cles), no major changes in the CC's second derivative is 
observed. 

In Fig. ^(b) the mean potential energy per particle, 
V, as a function of the total system energy is displayed. 
For the lowest densities two different regimes can be rec- 
ognized. First, a steep increasing behavior for energies 
E < 0.5eo, that can be related to an increase of the 
mean interparticle distance and the appearance of inter- 
nal surfaces, as the number of atractive bonds decreases. 
Second, a saturating behavior for larger energy values, 
E > 0.5eo- For more dense systems, no such feature can 
be observed. A smooth behavior is displayed by V. 

Another interesting feature can be noticed looking at 
the root mean square deviation of partial energies (ki- 
netic or potential), ay,if, shown in Fig. |](c). For the 
lowest densities (filled triangles and diamonds) an abrupt 
decrease is observed at E ~ 0.5eo, whereas for the high- 
est considered density, p h = 0.8er~ 3 (filled circles), no 
trace of such behavior can be reported. For future refer- 
ence, it is worth noting that the transition in the mono- 
tonic character of the curve takes place at a density value 
p ~ p c = 0.35eo(filled squares). 

To attain a better understanding of the behavior of 
the magnitudes displayed in Figs. 0(b), and (c), we have 
followed two strategies: i) to study particle-particle cor- 
relations, ii) to analyze fragment mass distributions ac- 
cording to the already presented fragment definitions. In 
what follows we present results obtained within the first 
approach, and postpone to the next subsection the anal- 
ysis of fragment properties. 

In Figj2| the distribution of interparticle distances, Tij, 
is shown. Two density values already considered in FigTy, 
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namely p h and p l , are shown in panels (a) and (b) re- 
spectively. 

For each density, three total energy values, lower 
than(solid line), equal to (dotted line), and larger than 
(dashed line) E — 0.5eo, were considered (see caption). 
In both panels, a vertical line was included indicating the 
interaction cut-off radius. 
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FIG. 2. Normalized interparticle distance distribution cal- 
culated for p l and p h systems are shown in panels (a) and 
(b) respectively. The following three total energy values, 
E = 0.0,0.5,1.0eo, were considered in panel (a), whereas 
E = —0.50, 0.5, l.Oeo for panel (b). They are shown with, 
solid, dotted, and dashed lines respectively. Note the differ- 
ent scales used for m in both panels. 

In panel (a) (very diluted case) a well defined structure 
of interacting particle pairs can be seen at low energies 
(continuous line). The displayed peaks signal the pres- 
ence of a rather large self-sustained drop, with a first, 
second and even third nearest neighbors structure. As 
energy increases, this structure fades away, and an im- 
portant reduction of the number of interacting particles 
can be noticed. This occurs as a consequence of the spa- 
cious volume available , that is big enough to accommo- 
date small non-overlapping particle aggregates. In this 
transition, well defined q-space clusters appear, surfaces 
are produced and V tends towards a residual value (see 
bellow). 

A completely different behavior is observed for the 
denser instance, p h , shown in Fig. ||(b). In this case, no 
changes in the statistical distribution of is observed 
as the energy is increased. (Note that the bell shaped 
distribution associated with the 'trivial' noninteracting 
pair counting is superimposed over the peak structure. 
However, the presence of a first, second, and even a third 
nearest neighbors can still be traced). The container im- 
poses a severe volume restriction on the system, and even 



for high total energy values, each particle is confined in 
the attractive concavity of partners potential, between 
the repulsive core of nearest neighbors. 

This last observation can complete a general picture, 
within which the behavior of the mean potential energy, 
V, for the high density cases can be easily understood. As 
more energy is added to the system no structural tran- 
sition is allowed by the constraining volume. Particles 
can not escape from neighbors most attractive potential 
range and a smooth increased in V, related to the aver- 
age time spent in the most negative potential areas, can 
be seen (Fig 0(b)). 

The decreasing of cry , observed for low densities (filled 
triangles, and diamonds in Fig. |l|(c)), can be associated 
to the loss of 'attractive bonds' between particle pairs 
that takes place when a non interacting light cluster 
regime dominates. As an increasing number of particles 
stop interacting with one another, the system dynamics 
gets 'less chaotic' (see for a dynamical characteriza- 
tion of this system), and a decrease of cry is induced. 

On the other hand, for high densities, the region with 
the strongest non-lineal nature becomes the most rele- 
vant part of the interaction potential. As a consequence 
of that, an increase of the fluctuations in potential en- 
ergy between successive configurations as total energy is 
added can be expected (see the behavior displayed by cry 
in Fig 0(c)). 

It is worth noting that the presented picture can also be 
used to interpret recent results regarding the behavior of 
the maximum Lyapunov exponent, MLE, in constrained 
systems (see j28|). Moreover, a striking similitude be- 
tween the behavior of <jy(E,p) and MLE(E, p) can be 
noticed, comparing FigQc), and Fig.5 of Ref. g§. 



B. Fragments Inside the Volume 



In section [II, three fragment recognition algorithms 
were described. Each one of them make use of different 
correlation information in order to link particles into clus- 
ters, and then give different physical information about 
the system under study. 
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mass 

FIG. 3. MST cluster mass distributions calculated for den- 
sities: p = 0.01, 0.03, and O.lOer -3 are displayed in panels (a), 
(b), and (c) respectively. Three total energy values were con- 
sidered for each density, Etot = 0.0,0.5, and l.Oeo.They are 
displayed as full, dotted and dashed lines respectively 

In Fig. H the results of the MST algorithm analysis 
are shown. MST spectra were calculated for different 
energies (E — 0.0, 0.5, and l.Oeo, displayed as full, dotted 
and dashed lines respectively) for three system densities: 
p = 0.01,0.03, and O.IOct -3 . They are shown in panels 
(a), (b), and (c), of Fig.|| respectively. 

At low densities, panel ||(a), the system evolves from 
a heavy cluster dominated behavior at low energies, to- 
wards a light cluster dominated behavior, at high ener- 
gies. This reflects the fact that the constraining volume 
is large enough to allow the system to fragment in well 
defined drops as energy is increased. It is interesting to 
notice that results displayed in Fig. ||(a) correspond to 
full triangle symbols in Fig.^J(a), i.e. the one for which the 
CC displays a loop. It was argued that big fluctuations 
in kinetic energy should be expected. This is indeed the 
case and it is related to the fact that well defined surfaces 
appear in the system (see (§)). 

On the other hand, as density increases, the mass dis- 
tributions converge to a u-shaped one. Almost no spa- 
tially well separated structures (i.e. MST clusters) can 
be identified, aside from the trivial huge cluster that com- 



prises almost all the system particles. 

In Fig. [|, the effect of taking into account momen- 
tum correlation in the definition of clusters is displayed. 
There, a cluster analysis for a system of N — 147 par- 
ticles with p = 0.2a- 3 using ECRA, MSTE, and MST 
prescriptions is presented. The mass distribution func- 
tion, and the fragment internal energy as a function of the 
cluster mass, Ei nt (s), are shown in the first and second 
columns, respectively. Three total energy values were 
considered: E to t = —0.6 (panels (a)-(b)), 0.0 (panels(c)- 
(d)), and 0.5eo (panels (e)-(f)). Empty triangles, solid 
squares, and empty circles correspond to MST, MSTE, 
and ECRA results respectively. 




mass mass 
FIG. 4. Fragment mass distribution, and cluster internal 

energy values as a function of the cluster mass, are shown in 

the first and second columns respectively. Three system total 

energies are considered: E to t ~ — 0.6eo, panels (a-b), O.Oeo, 

panels (c-d), and 0.5eo, panels (e-f). Empty circles, filled 

squares, and empty triangles, correspond to ECRA, MSTE, 

and MST data, respectively. 

It can be seen that, for this density, the MST pre- 
scription finds just a big drop. On the other hand, both, 
the MSTE and ECRA algorithms, find non-trivial cluster 
distributions that show the expected transition from a u- 
shaped, towards and exponential behavior, as the system 
energy is increased. As a general rule, ECRA algorithm 
produces more bound clusters than the MSTE one. This 
is consistent with the claim that the MBDF are found by 
the ECRA algorithm. Due to this feature, from now on, 
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just the properties of the system according to its ECRA 
cluster structure will be considered. 



C. Transition Signals 

As seen in Fig.^J, in which a system at a p = 0.2a~ 3 
density was considered, the ECRA cluster distributions 
undergo a transition from a U-shaped spectrum towards 
an exponentially decaying one. Moreover, the one corre- 
sponding to the intermediate energy value |^(b), displays 
a power law like behavior. In fact, for any other density 
studied, the same behavior can be reported, and a total 
energy value can be determined, for which a power-law 
like mass distribution can be found. 

This result is interesting, because it implies that some 
kind of transition between two regimes, one with, and 
other without large ECRA drops, i.e. with or without the 
presence of liquid-like structures, is taking place in high 
density systems. Moreover, this transition can not be 
detected simply using spatial correlations information. In 
this way, considering the appropriate cluster prescription, 
power-law mass distributions can be used to trace the 
transition line in a (T, p) diagram. 

Within the Fisher Droplet Model [p4j| , when dealing 
with infinite systems, one expects to find such a scale- 
free distribution at a critical point, where a continuous 
transition takes place, and scaling hypothesis can be ap- 
plied. However, it has been reported ^,^,^l[ that sev- 
eral critical signatures also appear when first order phase 
transition are analyzed in small systems. In particular, as 
stated in Ref , in small systems the largest cluster gets 
a size comparable to the vapor fraction before dissapear- 
ing when the system crosses the coexistence line. There- 
fore there is an energy for which a pseudo-invariance of 
scale and a resemblance with critical behavior can be 
expected for small systems undergoing phase transitions 
that, in the thermodynamical limit, would be univocally 
classified as first order ones. 

In order to search for a power-law behavior in fragment 
distributions the following single parameter function was 
employed to fit fragment mass spectra. The contribution 
of the largest fragment was disregarded pq] : 
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n(A) = q A T ,with q (r) 



(7) 



n(A) is the number of fragments with mass number A, 
and qo is a normalization constant. 

The quality of the fitting procedure was quantified us- 
ing the standard \ 2 coefficient (see ||^] for details), and 
an energy value, E*, was associated to the best fitted 
spectra. In Fig.|](a) a typical \ 2 calculation is shown for 
ECRA clusters in a p — 0.10<j~ 3 system. From the figure, 
a value of E„ = 0.10 ± 0.05eo, can be reported. 
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FIG. 5. Transition signatures calculated for & p — O.lOcr -3 
system. The x 2 fitting coefficient, the second moment of the 
mass distribution, M%, and the normalized mean variance of 
the largest fragment, NVM, are displayed in panels (a), (b), 
and (c) respectively. 

Another useful observable to search for scale-free dis- 
tribution functions is the second moment of the cluster 
distribution M^ = -A 2 n(.A). As m the determination 
of t, the largest cluster is excluded from the primed sum. 
Mi is proportional to the compressibility kt, that, in the 
thermodynamical limit, presents a power-law singularity 
at a second order transition. For finite system the di- 
vergence is replaced by a maximum. In panel |5|(b) it is 
shown the estimation of M2(E), for a p — 0.10cr~ 3 sys- 
tem. It can be noticed that the maximum is located at 
the energy value for which the spectra is best fitted by a 
power-law like dependence. 

This kind of agreement is also achieved when the nor- 
malized mean variance, NVM, of the largest fragment 
mass, A, 



is analyzed. NVM is defined as: 



NVM 



< A 2 > 

max 



^ A max > 
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and it proved to be a robust tool for the characteriza- 
tion of transition phenomena in which an enhancement of 
fluctuations occurs j33|. As can be seen from panel ||(c), 
this signal, although slightly shifted within the working 
resolution, also peaks in the same region of the previous 
ones. 

The signal agreement reported for the p = 0.10c -3 
case in Fig.g, is also achieved for every density analyzed 
in this paper. This means that, in the whole density 
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range, power-law mass distributions for phase space de- 
fined ECRA clusters are found whenever large fluctua- 
tions take place in the system. 

In order to properly characterize the state of the sys- 
tem which displays power-law mass distributions, the val- 
ues attained by the r exponents must be analyzed, as a 
function of the density. This dependence is shown in 
Fig.^(a). For a true critical phenomenon, taking place 
in three dimensional systems, 2 < t < 3 is expected 
This condition is not satisfied for systems with 
p < 0.05c~ 3 . For this highly diluted systems, the ob- 
served free-scale distribution is not expected to survive 
the thermodynamic limit. It is not related to any con- 
tinuous transition, but arises as a finite-size effect. It is 
worth noting that p ~ 0.05<r -3 is the maximum density 
for which the respective caloric curve displays a loop, and 
then, negative c v (see Fig.|l|). 
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FIG. 6. Estimated values from ECRA cluster distributions 
for the critical exponents r and 7 are shown in panels (a), 
and (b) respectively. The expected values for the 3D-Ising 
universality class are displayed by a dashed line. 

On the other hand, for p > 0.05, the calculated t ex- 
ponents show a rather good agreement with the t = 2.21 
value (dashed line), expected for liquid-gas transitions. 

In a second order phase transition, the behavior of M% 
near a critical point can be described in terms of the 
critical exponent 7 pffl: 



M 2 (e) oc 



(9) 



where, e = (Ec — E)/E c , measures the distance to the 
critical point. This relation is valid for an infinite system 
in the limit e — > 0. As already mentioned, in a finite sys- 
tem Mi displays a maximum instead of the divergence 



predicted in equation (Q). Having this in mind, a cal- 
culation procedure introduced in Ref. J27J (7- matching), 
that takes care of finite size effects, was used to calculate 
the 7 exponent value for our system (see also Ref. |[|] 
for details). 

The results of the 7 exponent estimation is shown in 
Fig.^(b). The obtained values tends toward "fiiq- ga s = 
1.23 (dashed line) value, that is expected for a liquid- 
gas transition. It is worth noting that this convergence 
is achieved for densities p > 0.35eo (this value equals 
p c in Fig.|l|). For this density, a change in the behavior 
of the potential energy fluctuations as a function of the 
total energy was repor ted in Figj^(c). This can be as- 
sociated (see Sec. IV A) to the onset of the invariance of 
the statistical interparticle distance distribution (FigJ|), 
that takes place as a consequence of the imposed volume 
restriction, and precludes the occurrence of a continuous 
transition. 



V. PHASE DIAGRAM 

In the previous section several signatures of a change in 
the properties of ECRA fragment mass distributions were 
analyzed. At any given density, a system energy, E*(p), 
can be determined for which large fluctuations appear in 
the system, and a scale-free like mass distribution can be 
found. 
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FIG. 7. Density dependence of the system energy, E t , at 
which transition signals are detected. The dashed line is in- 
cluded to guide the eye. 

In Fig.fi], the dependence of E* with the system den- 
sity can be seen. It can be noticed that for p > 0.35a -3 , 
a rather constant value is attained for the transition en- 
ergy. A similar result was reported in Ref. Jll| , using the 
already presented MSTE cluster definition. 

Some insight about the properties of the system along 
the line depicted in Fig.^, can be gained from the analy- 
sis of the interplay between the mean internal potential 
energy per particle stored in ECRA clusters, Vi nt , and 
the mean inter-cluster interaction energy, Vi C . 
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Vi 



i,iec fc 



(10) 

(11) 



Vi„t stands for the mean potential felt by a particle in 
a cluster due to its interaction with the other members 
of the same cluster. On the other hand, V lc represents 
the mean interaction potential felt by a particle due to 
its interactions with particles belonging to the other clus- 
ters. This magnitudes are displayed in FigJ| for densities 
p = 0.01,0.20,0.50, and 0.80cr~ 3 in panels (a) to (d) re- 
spectively. 
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FIG. 8. Vint (empty circles) and Vt c (empty squares) as 
a function of the system total energy are shown. Densities: 
p = 0.01,0.20,0.50, and 0.80a -3 are displayed in panels (a), 
(b), (c), and (d), respectively 

The general trend displayed by Vint and Vi C is easy 
to understand. At any given density, for low energies, a 
big ECRA cluster can be found having a large binding 
energy. In addition, as no many other clusters beside the 
biggest one exist, V~i c ~ 0. At high energies, the ECRA 



partitions turn out to present a rather high multiplicity 
of light clusters. Consequently Vi nt — ► in this limit and 
Vi C absolute value increases. 

After a close inspection of Fig.pl one can realize that, 
for any given density, the energy at which Vi n t — Vi C , hap- 
pens to be E Sf {p) 1 i.e. the energy at which power-law like 
mass distributions appear. This means that scale-free 
ECRA mass distributions have the following property: 
a balance is established between the mean potential en- 
ergy a particle, that belongs to a given cluster, feels due 
to its interaction with other members of the same clus- 
ter and the one associated with its interaction with the 
rest of the particles in the system. This potential energy 
balance, that does not allow to distinguish contributions 
from inside or outside of a cluster, is reminiscent of the 
vanishment of the chemical potential and surface tension 
terms that takes place at the critical point within the 
Fisher droplet model framework pifl . 

Gathering all the information obtained so far, the 
phase diagram for our 147-particle Lennard- Jones drop 
can be built. In Fig.^ the resulting phase diagram is 
presented. The empty circles are (T*,p) points obtained 
from the aforementioned cluster analysis. The full line is 
an estimation of the coexistence line, as obtained from 
the analysis of the system specific heat, c„ (see Eq||). In 
this case, for densities for which the respective CC dis- 
plays a loop, T(E*) has been taken as the average tem- 
perature between the two values corresponding to the 
location of the c v singularities. For denser cases, the 
location of the maximum of the c v has been taken. No 
reliable estimation of the energy that maximizes c v could 
be achieved for p > 0.35er~ 3 . 
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FIG. 9. Phase diagram for a 147 particle Lennard Jones 
drop. Cluster-based and c„-based determination of the coex- 
istence line is marked with empty circles, and full line respec- 
tively. The dashed line is included to guide the eye. 

From the figure, it can be seen that the symmetry be- 
tween low and high densities, characteristic for the infi- 
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nite size limit, is lost in the finite size system. An almost 
linear raising branch, for p > p c , can be seen instead. For 
every density value, below the depicted coexistence line, 
a liquid-like ECRA structure can be identified, whereas 
a vapor behavior can be observed above it. 

It is interesting to notice that three density regions 
can be identified. The region labelled A (p < 0.05cr -3 ) 
presents the signals expected for a first order phase tran- 
sition occurring in a finite system: the corresponding 
caloric curve shows a loop that can be associated to a 
negative specific heat, and a structural transition can be 
recognized (Fig.||(a)). These region is the only one in 
which the available volume is large enough for the sys- 
tem to fragment into a set of non-overlapping drops. 

On the other hand, in region C (p > 0.35er -3 ), where 
the system is rather compressed, and MST algorithm rec- 
ognizes just one big fragment, a second order transition 
seems to occur in phase space, whenever the transition 
line is crossed. No anomaly in the respective caloric curve 
is observed, no qualitative changes in the configuration 
statistical properties is reported (FigJ|(b)), and the cal- 
culated critical exponents, r, and 7, are in good agree- 
ment with 3D-Ising (liquid-gas) universality class. 

In between those two regions, region B (0.05cr~ 3 < p < 
0.35fi~ 3 ) can be identified. The finite size of the system 
plays a major role in this density range. Even a sensi- 
ble t value suggests that physically meaningfull scaling 
properties are present in the system, the corresponding 
7 exponent values are too low to classify the transition 
as a continuous one. 



nents, r and 7, related to the displayed power-law distri- 
bution, a further classification in density ranges can be 
obtained. 

In the abovcdefincd region A, the obtained t values are 
too big to be related to a true critical behavior. For den- 
sities 0.05er~ 3 < p < 0.35cr~ 3 , eventhough the values of 
t exponents are in good agreement with the correspond- 
ing 3D-Ising universality class, the value attained by 7 
critical exponents came out to be too low. Finally, for 
densities above p ~ 0.35er~ 3 , both, r and 7, are quite 
close to the accepted values for the 3D-Ising universality 
class. 

This is of particular interest for the statistical model 
approach because by definition, at freeze out, an ensem- 
ble of well defined fragments is to be constructed. Such 
assumption implicitly locates the system under study in 
region A, where first order phase transition is to be ex- 
pected. It is then clear that in the frame of such ap- 
proaches regions B and C are excluded from the analysis. 

A natural continuation of this work is to analyze the re- 
lation between fragmentation inside the container volume 
and the corresponding asymptotic mass spectra when 
walls are removed. This work is currently under progress. 

Finally, we would like to mention that the behavior 
of the equation of state for small systems has been un- 
dertaken in Ref. ju| . The results obtained in that work 
differs from ours (critical exponents) but they were based 
on properties of MSTE fragments, which we have disre- 
garded favoring the ECRA approach because of the rea- 
sons stated right after FigJ|. 



VI. CONCLUSIONS 



In this paper we have undertaken the analysis of ther- 
modynamical properties of a finite classical system con- 
fined in a volume. Apart from the basic interest on such 
a problem, it might be relevant on the frame of the de- 
scription of fragmenting system according to statistical 
models. 

We have been able to find, that from a coordinate space 
stand point the equation of state of such a system is 
quite simple. There is a maximum value of the density 
(p ~ 0.05er~ 3 ) up to which the system can undergo a 
first order phase transition. In this region, the associated 
CC displays a loop and the thermal response function 
attains negative values. This feature allowed us to define 
a transition curve. For higher densities, there is simply no 
room enough to allow the system to develop well defined 
internal surfaces, i.e. only one big drop can be detected. 

When one turns into a description of fragments in 
which correlations in p-space are taken into account, a 
much reacher structure appears. Now, even for densities 
bigger than p ~ 0.05a -3 transitions from u-shaped mass 
spectra to exponentially decaying ones are signaled by 
the appearance of scale-free mass distribution of ECRA 
clusters. However, from the analysis of critical expo- 
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